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We report a combined theoretical and experimental study of electrical transport in weakly-coupled 
doped superlattices. Our calculations exhibit negative differential conductivity at sufficiently high 
electric fields for all dopings. In low-doped samples the presence of impurity bands modifies the 
current-voltage characteristics substantially and we find two different current peaks whose relative 
height is changing with the electron temperature. These findings can explain the observation of 
different peaks in the current-voltage characteristics with and without external THz irradiation in 
low-doped samples. From our microscopic transport model we obtain quantitative agreement with 
the experimental current-voltage characteristics without using any fitting parameters. Both our 
experimental data and our theory show that absolute negative conductance persists over a wide 
range of frequencies of the free-electron laser source. 
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I. INTRODUCTION 



Perpendicular charge transport in biased superlattices 
is dominated by resonances due to the alignment of en- 
ergy levels in different wells. These resonances yield dis- 
tinct peaks in the current-voltage characteristics as- 
sociated with negative differential conductivity (NDC) 
at fields above the peak. The instability associated with 
NDC causes the formation of electric field domains ||] 
as well as self-sustained oscillations in such structures 
[Q. While for strongly-coupled superlattices the elec- 
tronic minibands dominate the electrical transport 
in weakly-coupled superlattices the transport is due to 
sequential tunneling from one well to the next. (For a dis- 
cussion of the appropriate regimes see Refs. 0,0].) This 
situation has been already regarded in Ref. B for tun- 
neling between the lowest level and excited levels in the 
adjacent well. There the current is driven by the different 
occupation of the two levels and a maximum of the cur- 
rent occurs when the different level are aligned. Tunnel- 
ing between equivalent levels at low fields is slightly more 
complicated, as alignment occurs at zero field, where, of 
course, the current vanishes. The key point is the treat- 
ment of broadening of the states due to scattering which 
essentially determines the transport. This idea has been 
exploited to determine scattering rates by studying the 
transport between two quantum wells |9|-pd]|. In the ex- 
periments |l0| , |rT| impurity scattering was diminished by 
the use of remote doping which enabled to study electron- 
electron scattering rates. In contrast to this we focus on 
doped superlattices in the present paper. There impurity 
scattering at the ionized donors is an important scatter- 
ing process whose impact we will examine in the follow- 
ing. In a previous study [O a heavily doped sample was 



investigated and good agreement with experimental data 
was found. Here we perform a systematic study of the 
low-field transport in such structures for different doping 
densities. We find that the formation of impurity bands 
]l3| , p^ for low-doped samples causes a strong tempera- 
ture dependence of the current-field relation which may 
display a double-peak structure at low fields. 

If the superlattice is subjected to an external mi- 
crowave field, photon-assisted tunneling (PAT) is possi- 
ble where replica of the resonances are observed at biases 
which differ from alignment conditions |p5|-pT[ by integer 
multiples of the photon energy. For certain field strengths 
of the irradiation field absolute negative conductance has 
been observed experimentally p8| , |l9| . The main features 
of these experiments could be described qualitatively 
]l8|-p0| within the standard theory of photon-assisted 
tunneling [^,^ but modifications due to photon side- 
bands from a single quantum well p3t have also been 
suggested to explain the experimental findings. Here we 
present additional experimental data and show that full 
quantitative agreement between theory and experiment 
can be found by a combination of a microscopic trans- 
port model with the standard theory of photon-assisted 
tunneling. This comparison strongly supports our claim 
that a microscopic treatment of impurity scattering is 
necessary for a full understanding of transport in low- 
doped superlattices. 

The paper is organized as follows: Our transport model 
is presented in Section II. In order to understand the 
generic behavior we give a phenomenological approxi- 
mation where many features can be seen analytically in 
Section III. The calculated results for different doping 
densities are presented in Sections IV and V using dif- 
ferent screening models, respectively. Our calculations 
are compared with two different experiments concerning 
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a highly-doped and a low-doped sample in section VI. In 
section VII we consider transport under external irradi- 
ation. Finally, we will discuss the general significance of 
our results. 



II. THE MODEL 

We consider weakly-coupled semiconductor quantum 
wells of period d. Then the electrons are essentially lo- 
calized in the wells and a reasonable basis set of wave 
functions is given by a product of Wannier functions 
'i'''{z — jd), which are maximally localized in well 
J, and plane waves e*'^ "". The z direction is defined to be 
the growth direction and k, r are vectors within the (a;, y) 
plane, ly denotes the subband within the well. Here we 
restrict ourselves to the lowest level and omit the index 
ly in the following and the energy of the lowest level is 
used as a reference point. 

Regarding only next-neighbor coupling Ti we have the 
following Hamiltonian {F is the electric field, and e < 
is the charge of the electron) : 



H = 



Y,[iEk-3eFd)a]{k)a,(k) 



Tiflt (k)a,(k)+Tia^(k)a,+i(k) 



(1) 



with the in-plane kinetic energy Ek = ?i^A;^/(2m.u,), 
where m^, is the effective mass in the well, a., and 
are the annihilation and creation operators of electrons 
in well j, respectively. denotes the contribution 

due to scattering which is not k-conserving. 

Within the lowest order in the coupling Ti the current 
density from the lowest level in well j to the lowest level 
in well j + 1 is given by 

2e,^, ,^ f°° dE , 
J.^..^ = ^T.\^.fJ_^^A,,.iKE + eFd) 

xAj{'k,E)[nFiE - ^j) -npiE + eFd- fij+i)] . (2) 

Here e is the electron charge, A is the sample area, 
fij is the local chemical potential in well j measured 
with respect to the energy of the lowest level. npiS) = 
-I- exp(f /fcsTe)], and Tg is the electron temperature. 
Fd denotes the voltage drop per period d. The spec- 
tral function >l(k, E) is calculated for a given intrawell 
scattering H^'^'^^^ via the retarded self-energy E™*(k, iJ): 



A{-k,E) = 



-2Im{I]'°'(k, E)} 



[E -Ek- Re{E'°'})^ + (Im{S™t})' 



(3) 



/ij is related to the electron density Uj in well j via the 
relation 



/oo 
dE pj {E)nF{E — jij) 



with the density of states 
2 



2t:A 



(4) 



(5) 



where the factor 2 reflects the spin degeneracy. 

In our microscopic calculation we proceed as follows: 
First we determine the coupling Ti as well as the Wan- 
nier functions ^(z — jd) for the given superlattice pa- 
rameters (see Appendix ^). Then we calculate the self- 
energy (k, E) for impurity scattering using the self- 
consistent single-site approximation shown in Fig. |l|. The 
respective formulas are given in appendix The ma- 
trix element for impurity scattering is calculated from 
the Coulomb potential of the individual ionized donors. 
Screening is treated in two different approaches, the 
random-phase approximation (RPA) for a free-electron 
gas and the Thomas-Fermi approximation (TF) using the 
actual density of states at the Fermi level (see appendix 
^). Using the calculated spectral functions A{k,E) the 
chemical potential is determined by setting rij = Njj in 
Eq. (^), where No is the doping density per period. Fi- 
nally the current is calculated from Eq. (^. Note that 
all quantities used in the calculation are defined by the 
sample parameters and no fitting parameters are used. 



III. PHENOMENOLOGICAL DESCRIPTION 

In this section we want to provide some insight into 
the question how scattering effects the transport. Using 
a constant self-energy we derive some simple expressions 
for the current-field relation which will help to under- 
stand the full calculations presented in subsequent sec- 
tions. 

As mentioned in the introduction the level broaden- 
ing essentially determines the transport in the sequen- 
tial limit. This can be easily seen in the limit of van- 
ishing scattering. Then the spectral functions become 
(5-functions, A(k,E) = 2ttS{E - Ek). In this case the 
current vanishes for eFd ^ 0. (In addition further reso- 
nances may occur at finite fields, when the lowest level is 
aligned with higher levels in the neighboring well, which 
are not considered here). Assuming eFd > we rewrite 
Eq. (^) as follows 

J,^,+i j^^dE{A,+,)iE,F)p,{E) 

X [npiE - p.j) - np{E + eFd - pj+i)] (6) 

with 

Jo^dE,A,{KE)A,+,ik,E + eFd) 



J,dEkA,{k,E) 
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where we used Eq. ^ and performed the contin- 
uum hmit. Now let us assume a constant self-energy 
E''<'*(k,£;) = -iT/2 in Eq. (0) for the sake of simplic- 
ity. Then the spectral functions become Lorentzians 
A{k,E) = r/[{E - Ekf + rV4]. Extending the lower 
integration to — cx) we obtain 



{A. 



2T 



j+ii 



{eFdf + r2 



(8) 



which only depends on F. Note that this simple model 
with a constant self-energy cannot be used in the calcu- 
lation of the density of states (||) as the integral for the 
electron density (^) diverges in this case. Therefore we 
use the free-electron density of states Pj{E) = poQ{E) 
(with pq = m/iTrfi^)) and obtain for equal chemical po- 
tentials — fJ.j+1 — fJ.: 



J{F) = epo 



2T 



h (ei^d)2+r2 



'.Fd 



dEnpiE-n). (9) 



For low electron temperature and voltage drop 
{ksTe, eFd < p.) we find 



J{F) = epo 



2reFd 
T (eFd)2 + r2 



(10) 



Thus, we recover an ohmic behavior for low fields eFd <^ 
r, a maximum of J{F) at eFd = T, and negative dif- 
ferential conductivity for eFd > T. Eq. ( [lo| ) has been 
essentially used in Refs. ^0 11 1 for the determination of 
scattering rates T/H from tunneling between two two- 
dimensional electron gases. Similar models using a phe- 
nomenological broadening T have been applied to the se- 
quential tunneling in superlattices in Refs. po| , ^ . The 
current at the maximum is given by 



J 



epo 



(11) 



which is independent of doping, scattering, and temper- 
ature in the limit of p » eFd, ksTe considered here. 

If fcsTe becomes of the order of p the factor 
J^^"^ dE npiE — p) in Eq. (||) is smaller than eFd and 
we obtain a decrease of the current with temperature. 
Here we have to take into account the temperature de- 
pendence of the chemical potential p. From Eq. (^) we 
find 1 -I- exp(/x/fcsT'e) = eicp{n/ poksTe). This gives a 
zero field conductivity 



dJ _ 2e^po dTl 
dF\F=o ~ T~ 



exp 



PoksTe 



(12) 



which is almost constant for ksTe < n/po and drops as 
l/ksTe for large temperatures as observed experimen- 
tally in Ref. [Q. For completeness, we give the result in 
the high temperature limit (fc^Tg 3> eFd,n/po) 



J(F) = 



en 



2VeFd 



keTe n {eFdf+V^ 



(13) 



which follows directly from Eq. (g). It is interesting to 
note, that Eq. ( |l0[ ) is identical to the current-field re- 
lation calculated for miniband conduction using a 
constant scattering time Ti/T for p > 2Ti and ksTf, — 0. 
Equivalently, Eq. ( p^ has been obtained from miniband 
transport in the limit ksTe ^ 2Ti,n/po as well 
This shows that the models of sequential tunneling and 
miniband conduction give the same results provided ei- 
ther the electron temperature or the electron density are 
large. 



IV. RESULTS FOR RPA SCREENING 

As a model system we choose an Alo.sGao.yAs-GaAs 
superlattice with barrier width h = 10 nm and well 
width w = \Q nm. We use the conduction band offset 
240 meV and the effective masses ruw = 0.067me and 
TOfc = 0.0919me |3^] in the Kronig-Penney model yield- 
ing a couphng Ti = —0.0116 meV. We assume J-doping 
in the middle of the quantum wells. The interaction with 
impurities located in different wells is found to be neg- 
ligibly small. In this section screening is treated within 
the RPA assuming a free-electron gas. 



A. Density of states 

In Fig. H we show the resulting densities of states for 
four different doping densities Nd . For high Nd the den- 
sity of states exhibits a monotonic increase from p = at 
E < -Emin to p ^ po for E ^ (X, where i?min denotes the 
lowest edge of the density of states. In contrast to this the 
density of states splits into two parts for small doping: 
p{E) takes finite values in a certain region below E — 0, 
which we will refer to as an impurity band. For higher en- 
ergies p{E) is quite similar to the density of states of the 
free-electron gas. These results are in good agreement 
with the findings of Ref. |l^ . The onset of the impurity 
band occurs at slightly larger energies i?min here, as the 
wavefunctions are less confined due to the spreading into 
the barrier which was neglected in Ref. . 

We also marked the positions of the Fermi level E^ 
(i.e., the chemical potential for — 0). For low densities 
the position is just in the middle of the impurity band, 
indicating that the impurity band consists of exactly two 
states per impurity due to the assumed spin degeneracy. 
(This degeneracy would be lifted if spin-resolved interac- 
tion was taken in account, see also appendix^.) For high 
densities the position of E^ roughly equals the Fermi 
level of the free-electron gas No/po- The crossover be- 
tween these two limits occurs at Nd sa S-lO^^/cm^ where 
E^^^O. 

The respective spectral functions are plotted in Fig. ||. 
For E — 5 meV A{k, E) resembles a Lorentzian centered 
close to Eh ~ E. This is the generic behavior of a free 
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quasiparticle with a finite lifetime due to scattering. The 
width of the spectral functions is increasing with doping 
due to the enhanced scattering. We find a full width at 
half maximum F ~ 0.5 meV for Nd « 1 • 10^'^/cm^ and 
r = 5 meV for « 1 • 10^^/cm^, which arc in the range 
of the calculated values of — 2Im{I](k, E)}. 

For E = —5 meV the spectral functions exhibit a 
monotonic decrease. For high doping the slope is com- 
parable to the slope at i? = 5 meV. In contrast to this, 
the spectral function for i? = 5 meV and E = —5 meV 
are entirely different for low doping, indicating that two 
different types of states occur. While the states for 
= 10^°/cm^ are essentially free-particle states at 
E ^ 5 meV, they are localized in space for = — 5 
meV, which is the signature of an impurity band 



states for fc^Tg > li^minl (see Fig. ||(b) where the grey 
scale denotes the relative occupation). The current con- 
tribution due to the free-electron states can be under- 
stood within the phenomenological constant-F approach. 
There is a maximum at eF\o^d « F, which coincides with 
the full width at half maximum of the spectral function 
aX E — b meV in Fig. ^. The amplitude of this maximum 
depends on two competing effects: On the one hand the 
occupation of the free-electron states increases with tem- 
perature. On the other hand the Fermi-factor in Eq. (||) 
strongly decreases with temperature. This explains the 
calculated behavior, where the peak at eF\o^d takes its 
maximum at intermediate temperatures. 



V. RESULTS FOR THOMAS-FERMI SCREENING 



B. Currents 

We calculate the current densities Jj^j+i {eFd) for dif- 
ferent electron temperatures T^, from Eq. (^). The results 
are shown in Fig. ^. For all temperatures and densi- 
ties we find an ohmic range for low electric fields and 
negative differential conductance for high electric fields. 
Let us first regard the high doping case (a-c), where 
no impurity bands form and the Fermi level is signif- 
icantly above E — 0. In this case the approximation 
is justified and indeed we find a maximum at val- 
ues of eFd which are in the range of calculated values 
of F = -2Im{I]'°*(k,£')}. The height is estimated by 
^max = 0.91A/cm^ from Eq. ( [n] ) which is in good agree- 
ment with the full calculation at Tg = 4 K. Note that 
the maximal current is almost independent on the dop- 
ing in this range. For Nu = 5 • 10^^/cm^ the chemical 
potential is larger than ksTe for all temperatures. Thus, 
the current is hardly affected by the temperature. In 
contrast to this the current drops with temperature for 
lower doping {Nd ~ 10^^/cm^). All these findings are in 
good agreement with the phenomenological description 
using a constant F discussed above. 

For the low-doped samples (see Figs. ^(e,f)) an entirely 
new scenario occurs. Here we find two different max- 
ima in the current-field relation whose relative weight is 
changed by temperature. The reason for this behavior 
is the presence of impurity bands for these doping lev- 
els. For Te = 4 K we find a maximum at ei^high^ ~ 8 
meV. This is due to tunneling from the impurity band 
to the free states (see Fig. ^(a)). The maximum occurs 
at the energy where the bottom of the impurity band in 
one well is aligned with the band edge of the free-electron 
states in the neighboring well, i.e., eFhighd ~ l^^minl- An 
increasing temperature leads to a transfer of electrons 
from the impurity band to the free-electron states and 
consequently the current at ei^highc' decreases with in- 
creasing Tg. The density of states in the impurity band 
is much lower than in the free-electron states, and hence 
the majority of the electrons will be in the free-electron 



The properties related with the formation of impurity 
bands are sensitive to the actual screening of the interac- 
tion For low doping densities the density of states 
differs essentially from the free-electron density of states 
and thus the use of RPA-screening by a free-electron gas 
is questionable. In order to take this effect into account 
we use the Thomas-Fermi approximation (TF) with the 
actual density of states at the Fermi level (see appendix 
^) in this section. Of course neither the free-electron 
RPA nor the TF approximation treat the screening en- 
tirely correctly, but we hope to obtain some insight into 
the general features by comparing these two approaches. 

In Fig. ^ we show the resulting density of states which 
is in qualitative agreement with the results of the RPA 
screening (Fig. I). For Nd^5- lO^/cm^ the density of 
states is almost identical while for lower densities some 
deviations occur. Especially the onset of the impurity 
band Smin is shifted to lower energies for TF-screening. 
Furthermore the impurity bands extend over a larger en- 
ergy range and have a lower density of states, so that the 
total density is conserved. The reason for these devia- 
tions lies in the fact that TF-screening is less effective 
than RPA screening if the actual density of states at the 
Fermi level is used. Therefore both the binding energy 
of the impurities as well as the broadening of the states 
become larger. 

This manifests itself in the calculated current densities 
(see Fig. For high doping (a) the characteristics are 
almost identical, while for lower doping deviations occur. 
At first note that the maxima due to tunneling between 
free-electron states (the maximum for Njj — 10^^/cm^ 
as well as the maxima Flow for Nb = 2 • 10^"/cm^ and 
Nu = 10^°/cm^) are shifted to the right according to 
the stronger scattering which increases F. At second the 
peak at Fhigh is shifted to the right compared to Fig. ||. 
Again we find eFhighrf ~ E^i^ for both densities. 

Therefore we conclude that within both approxima- 
tions for screening the two maxima are determined by 
specific quantities describing the scattering. eF\ov,d re- 
flects the average broadening F of the free-particle states 
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and C-Fhighrf is the energy separation between the onset 
of the impurity band and the free-particle states. 



VI. COMPARISON WITH EXPERIMENTS 



Previously |12(| , the formalism was applied for the 
highly-doped sample {No = 8.75 • lO^/cm^) of 
Refs. [^,|2|. Good quantitative agreement with the ex- 
perimental data was found, albeit using a barrier width 
being 10% smaller than the nominal value. (A similar 
width had been used in the original analysis by the ex- 
perimentalists as well [Q.) The position of the first max- 
imum occurred at eFd =13 meV, which is almost inde- 
pendent of the barrier width (which mainly changes Ti) 
and in excellent agreement with the experimental finding. 
The second resonance, as well as the formation of field 
domains was also studied in Ref. [|2j, and again good 
agreement with experimental data was found. 

A low-doped superlattice (A'^o = 6 x 10^/cm^, b — 5 
nm, w = 15 nm, A = 8/im^) with A^ = 10 wells was 
used in the experiments of Refs. p8| , p^ in order to study 
the transport under strong THz irradiation from a free- 
electron laser. Additional data for this sample will be 
given in the following. Without irradiation a broad max- 
imum was found in the range 50 mV< C/ < 100 mV 
where the current is almost constant. For U > 100 mV 
domain formation sets in. Dividing by the number of pe- 
riods {N = 10), the maximum extends to ei^unirrC^ ~ 10 
meV. In contrast to this, the photon-replica under strong 
THz-irradiation could be consistently explained by as- 
suming an "instantaneous" current-voltage characteris- 
tic with a distinct maximum at [/ « 20 mV (i.e., 
eFii-rd = 2 mcV). 

We calculate the current-field relation for this super- 
lattice using the experimental sample parameters. In or- 
der to model the homogeneous doping we use 8 equally 
spaced (5-doping layers per period. The calculated den- 
sity of states for both RPA and TP screening is shown 
in Fig. ||. The density of states resembles the result for 
low doping found before. Nevertheless we do not find 
a separation between the impurity band and the free- 
particle states. The reason is the homogeneous doping: 
The different impurity positions have different binding 
energies which smears the impurity band. Again the on- 
set of the impurity band occurs at significantly lower en- 
ergies within the reduced Thomas-Fermi screening. Both 
values of l£'min| are smaller than the corresponding val- 
ues for low doping for the calculation done before (see 
Figs. HJI) . This is due to the larger well width in the sam- 
ple: The Wannier-states are less localized and therefore 



the matrix-element for impurity scattering(B2) as well as 
the binding energy of the impurity levels is smaller. 

The results for the current-field relation are shown in 
Fig. Again we find two maxima whose relative height 
changes with temperature. The position of the maximum 
for low temperatures, ei^highd, is almost identical to the 



value of |i?min| for both types of screening like in the 
calculations shown before. 

Now we can offer an explanation for the two different 
maxima occurring in the experiment |^,^ with and 
without irradiation mentioned above. For low electron 
temperatures and without irradiation the maximum at 
e^iiighC^ dominates the transport and thus domain forma- 
tion sets in at voltages exceeding U « A^ei^highC^ where 
A^ = 10 is the number of wells. If the THz radiation is 
present the electrons are excited from the impurity band 
into the free-electron states corresponding to a larger 
effective electron temperature. Thus, the maximum at 
U — NeFio^d is dominant, and the photon replicas cor- 
responding to this feature are seen experimentally. The 
experimental values therefore suggest e-FhighC^ — 10 meV 
and eF\ovid = 2 meV which is in excellent agreement with 
the calculation using Thomas-Fermi screening. This indi- 
cates that that the RPA using a free electron gas overesti- 
mates the screening in low doped samples. The Thomas 
Fermi approximation with the actual density of states 
at the Fermi level seems to reproduce the experimental 
results better in agreement with our argumentation in 
appendix |^. Therefore we will use it in the following for 
comparison with the experiment. 

In Fig. |l^ we compare the calculated currents with 
the experimental current-voltage characteristic without 
irradiation in a wider range of fields. Here we included 
the resonance around eFd « 50 meV between the lowest 
level and the first excited level as well. The calculation 
of the corresponding current is completely analogously to 
Eq. (|), for details see Ref. (t). Note that there are no 
fitting parameters involved in the calculation - all quan- 
tities including matrix elements and spectral functions 
are directly calculated for the given sample parameters 
as outlined above. Let us first focus on the low field re- 
gion. For t/ < 10 mV the experimental data are in good 
agreement with the calculated currents for = 4 K, 
the experimental lattice temperature. With increasing 
bias, the experimental data leave the 4K curve and fol- 
low the Te = 35 K curve at the plateau between 50 mV 
and 100 mV. This can be understood by electron heating 
inside the sample: For a voltage drop of 8 mV per pe- 
riod and a current of 0.6^A the Joule heating is P w 10 
pW per electron. In a recent transport experiment a 
distribution function with Tg « 40 K was observed p^ ] 
for this amount of heating albeit using a sample with 
higher doping. This shows that the electron tempera- 
ture strongly deviates from the lattice temperature in 
the experiment considered in good agreement with our 
findings. At [/ > 100 mV, the homogeneous field distri- 
bution becomes unstable, as the region of negative differ- 
ential conductance is reached and electric field domains 
form, causing the saw-tooth shape of the characteristic 
(see Ref. and references cited therein). For U > 450 
mV one can clearly see the resonance between the lower 
level and the first excited level in the well which is located 
48 meV above the ground level. Again the calculation 
exhibits two different peaks depending on the electron 
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temperature due to the different occupation of the im- 
purity bands, although only the high temperature result 
should be meaningful due to the heating of the carriers. 
The peak height of around 14/iA is in excellent agree- 
ment with the value of 13.6/iA found experimentally for 
our sample. The experimental peak position is located 
at a higher bias. This may be due to a voltage drop in 
the receiving contact, where a low-doped spacer layer of 
^contact = 50 nm thickness exists. If the electric field in- 
side the sample is large, it cannot be screened within the 
spacer layer and the effective field inside the sample is 
U ^ [Nd + dcontact)F instead of [/ = NFd as assumed 
in the figure. 

In order to circumvent the problems of electron heating 
we have investigated the temperature dependence of the 
zero-bias conductance G = dl/dU, where should be 
equal to the lattice temperature T. The results are shown 
m Fig. |ll| both for our full calculation using TF-screening 
as well as for spectral functions calculated within the 
self-consistent Born-approximation (B7) where no impu- 
rity bands form. In the latter case G is monotonously 
decreasing in T as shown in Fig. |Tl|. This can be eas- 
ily understood within the phenomenological constant-F 
approach (p^. However, a different scenario emerges if 
the electrons occupy impurity bands for low tempera- 
tures. Then G is strongly suppressed due to the small 
values of A(k, E) for E < 0, see Fig. As tempera- 
ture is increased, more electrons are excited to the free- 
electron states, and G increases with T until the impurity 
bands are almost empty at fc^T ~ |£'min|- This physical 
picture is confirmed by the experimental data shown in 
Fig. At low temperatures the agreement is quanti- 
tative, while at intermediate T the theory overestimates 
G; this is most likely due to additional scattering pro- 
cesses not included in our calculation, or by the presence 
of a contact resistance which may limit the experimental 
conductance. 

Thus we may conclude that the results of our calcula- 
tions are in good agreement with experimental data both 
for high and low doping. Nevertheless, a direct observa- 
tion of the two-peak structure is not available so far. 



VII. PHOTON-ASSISTED TUNNELING 

The standard theory of PAT considers tunneling be- 
tween two reservoirs between which a dc- voltage AJ7dc is 
applied. Let us denote the current-voltage characteristic 
(IV) without irradiation by I'^'^{AUdc)- Under irradiation 
an additional ac-bias U = A/7ac cos(27ri^t) is induced be- 
tween the reservoirs. Then the time-averaged IV is given 
by [||: 

QO / Ih \ 

/'"(A[/dc) = iMo')?!"" ^Udc + — ) > (14) 

l=-oc ^ ^ 

where a = eAUg.c/{hi') and J; is the ordinary Bessel 
function of order I. Thus, the current under irradiation 



is given as a sum over the photon replicas eA/7dc + Ihv 
where the alignment of the wells is shifted by integer 
multiples of the photon energy. The great advantage 
of Eq. ( |l^ is that it expresses all transport proper- 
tics in terms of I'^'^{AUdc)- Eq. ( |l^ has been apphed 
to photon-assisted tunneling in weakly-coupled superlat- 
tices by identifying AUdc—Fdcd and AC/ac = ^acc', where 
-Fac is the field component of the microwave field in the 
growth direction of the superlattice . Note that 

Eq. (|lj) also holds within a miniband model for strongly- 
coupled superlattices p4) . 

In Refs. |l8|j2C| ] the I'^'^ (eFdcd) curve was calculated 
phenomenologically and a qualitative agreement with the 
experimental data could be obtained. To refine the the- 
ory we use the results of our microscopic calculation (see 
Fig. |l^) here. As the external irradiation heats the elec- 
tronic distribution for zero bias as well, we use the curve 
for Te = 35 K. At this electron temperature about 50% 
of the electrons are occupying the states in the impurity 
band. Of course the actual electron distribution may 
deviate from a Fermi distribution under the strong irra- 
diation. Nevertheless we expect that an effective temper- 
ature approach gives a reasonable description of the ex- 
cited carriers. Further calculations show that the results 
for Te = 30 K or Tg = 40 K do not differ qualitatively. 
Quantitative agreement between theory (Fig. |l2|a) and 
experiment (Fig. p^) is found for hi^ = 6.3 meV (1.5 
THz) for different strengths of the laser field. We find 
a direct tunneling peak at Udii — NF\o„d 20 mV and 
photon replicas at J7 « Udir + Nhv/e = 83 mV and 
U « Udh + 2Nhv/e = 146 mV. Photon-rephcas of the 
second peak around A^-Fhighrf ~ 100 mV are less pro- 
nounced as this peak is broader. Our calculations show 
that they become visible if larger photon energies are 
used. For low bias and high intensities there is a region 
of absolute negative conductance |lj] , which we focus on 
in the following. 

In Fig. |l^(d) the laser intensity has been tuned such 
that maximal absolute negative conductance occurred for 
each of the different laser frequencies. Then we observe 
a minimum in the current at C/ « —Udn + Nhv/e which 
is just the first photon replica of the direct tunneling 
peak on the negative bias side. This replica dominates 
the current if the direct tunneling channel is suppressed 
close to the zero of Jo(a) in Eq. (p^), i.e., a w 2.4, as 
used in the calculation (Fig. [T^c). Both the theoreti- 
cal and experimental results show that absolute nega- 
tive conductance persists in a wide range of frequencies 
but becomes less pronounced with decreasing photon en- 
ergy. In the calculation absolute negative conductance 
vanishes for hv < l.d> meV which is approximately equal 
to hv < eF\ov,d. The latter relation has been verified by 
calculations for different samples as well. For hi/ = 5.3 
meV a smaller value of a = 2.1 (thin line) agrees bet- 
ter with the experimental data (in the same sense the 
value a = 2.0 agrees better for hv = 6.3 meV, com- 
pare Fig. p^). This may be explained as follows: If 
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strong NDC is present in doped superlattices the homo- 
geneous field distribution becomes unstable and either 
self-sustained oscillations or stable field domains form 
p5| . Then the IV deviates from the relation for homoge- 
neous field distribution, where U — NFd, and typically 
shows less pronounced NDC. Therefore maximal negative 
conductance is observed at a laser field corresponding to 
a value of a < 2.4, where the NDC is weaker and the 
field distribution is still homogeneous. The presence of 
an inhomogeneous field distribution could explain the de- 
viations between theory and experiment for U > 150 mV 
as well. 



VIII. DISCUSSION 

We have investigated the electrical transport in 
weakly-coupled doped superlattices, where the transport 
is given by sequential tunneling. Our calculations give 
negative differential conductance for all doping densities 
and temperatures at sufficiently large electric fields. This 
will give rise to instabilities leading to domain formation 
p6| , p7[ or self-sustained current oscillations Q. Within 
the full transport model using Eq. (||) these effects are 
discussed in Ref. 0. 

For high doping Nd ^ lO^'^^/cm^ or high temperatures 
ksT ^ |i?min| the electrons mainly occupy free quasipar- 
ticle states. Then the general behavior can be under- 
stood within a phenomenological model using a constant 
self-energy —iT/2. The current exhibits a maximum at 
eFd « r which can be used to investigate scattering 
processes. For doped samples impurity scattering is an 
important scattering process which we considered here. 
The inclusion of further scattering processes like inter- 
face roughness scattering, electron-electron scattering, or 
phonon scattering will increase F and therefore the posi- 
tion of the first peak. 

For low-doped samples ^ 10^^/cm^ and low tem- 
peratures ksT ^ I-Emini the presence of impurity bands 
influences significantly the low- field transport. Then a 
second maximum at eFhighd w l-Eminl occurs. This max- 
imum provides a possibility to obtain information about 
the position of the impurity band. This position depends 
strongly on the screening as can be seen by comparison 
of the calculations within RPA and TF. Therefore such 
experiments could serve as a test on various models for 
screening. For the sample considered here i?niin calcu- 
lated within the free-particle RPA is too low compared 
with the experimental onset of domain formation. In con- 
trast the result using TF screening gives excellent agree- 
ment with the experimental data. This indicates that 
screening within the free-particle RPA is too strong for 
low-doped samples. 

Furthermore it would be interesting to see if effects 
due to spin splitting of the impurity band are visible in 
experiments. These experiments can be both carried out 
in doped superlattices as well as in resonant tunneling 



between neighboring two-dimensional electron gases in 
the spirit of Ref. ||l^,|ll[. The latter has the advantage 
that problems due to domain formation in the region of 
negative differential conductivity are absent. An impor- 
tant aspect in such experiments is the problem of electron 
heating as the temperature refers to the temperature 
of the electronic distribution. In order to avoid heating, 
structures with thick barriers should be used where the 
Joule heating becomes small. 

For external irradiation we have demonstrated both 
experimentally and theoretically that absolute negative 
differential conductance persists in a wide range of fre- 
quencies hv > eF\ovid. The calculated current-voltage 
characteristics are in excellent quantitative agreement 
with the experimental data using a microscopically calcu- 
lated I'^'^ieFd) combined with the Tucker formula (|lj). 
Recently, the same model has been applied to photon- 
assisted tunneling in a different sample fs^ and quanti- 
tative agreement has been achieved as well. This shows 
that the simple Tucker formula allows for a quantita- 
tive description of photon- assisted transport in weakly- 
coupled superlattices. 
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APPENDIX A: CALCULATION OF THE 
TRANSITION ELEMENTS 



In a superlattice structure the coupling between neigh- 
boring wells Ti is related to the dispersion relation E{q) 
of the miniband (see [Q) via 

Ti = A dqEiq)e^'^^ . (Al) 

For a next-neighbor tight-binding model we have E{q) = 
— {A/2) cos{qd) and Ti is equal to a fourth of the mini- 
band width A. Here we calculate E{q) for a given 
superlattice via the Kronig-Penney model for the respec- 
tive sample parameters. Furthermore we determine the 
Wannier-functions ^'(z — jd) localized in well j from the 
Bloch- functions (t)q{z), where we choose the phase of the 
Bloch functions such that the Wannicr functions are max- 
imally localized [^. These Wannier-functions are used 
for the calculation of the matrix elements for scattering. 
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APPENDIX B: CALCULATION OF THE 
SELF-ENERGIES 

We assume that the electron density in the conduction 
band is provided by doping of the superlattice. Thus, 
scattering at ionized impurities is an important scatter- 
ing process. In addition there may be interface roughness 
scattering, phonon scattering, or electron-electron scat- 
tering, which we will neglect in the following. For weakly- 
coupled superlattices, the dominating scattering process 
occurs within the wells, which are assumed to be identi- 
cal. Thus, the well index j can be omitted. Scattering at 
the ionized impurities is described by the Hamiltonian 



k,p,a 



(Bl) 



where the subscript a denotes the impurity located at 
the position (r^jZa). The matrix element is calculated 
with the Wannier functions yielding: 



47res eo y/\r - r^p + (z - Za 

o2 



2eseop 



dz \E'*(z)*(z)e 



-pIz-ZqI -ip-i 



(B2) 



Using the bare Coulomb interaction Vq(p) the relevant 
integrals in the self-energies are divergent. Thus screen- 
ing is essential for the calculation. We treat screening 
within the Random-Phase approximation (RPA) of the 
free-electron gas as well as within an effective Thomas 
Fermi approximation(TF) (see appendix |c|) . With the 
screened impurity potential V^'^{p) — Va{p)/e{p) the 
self-energy is calculated within the self-consistent single- 
site approximation (shown diagrammatically in Fig. |^) 
like in Ref. ||l3|| . Then the self-energy contribution from 
the impurity a is given by 

E„(k, E) = j^Y. ^""(k - ki)G(ki, i?)F-(ki - k) 



P E l^-(k-ki)G(ki,i?)F-(ki-k2) 



A3 



ki,k2 



X G(k2,^)^„^^(k2 -k) 



where Gik,E) = {E - Eu - 
tarded Green function and 



(B3) 

i;"*(k, is the fuU m- 



S''^*(k,i?)=EE„(k,i?) 



(B4) 



is the sum over all contributions. In case of 5-doping 
the impurity potentials (B2) from different impurities lo- 
cated in the same well differ only by a phase factor and 



the sum over a can be replaced by a multiplicat ion with 
the number of impurities per layer NdA. Eq. (B3) can 
be transformed to the self-consistent equation (see, e.g., 

ill) 

if«(ki,k,S) = -K:^(ki-k) (B5) 

+ i E ^"(''l - k2)G(k2, i?)i^„(k2, k, E) 

which we solve numerically for a given self-energy func- 
tion ^^'(k, i?) entering G(k2,i?). We parametrize ki,k 
by Ek^ , Sfe, and = Z(ki, k) and discretize the resulting 
equation. This gives a set of linear equations for the com- 
ponents of K{Ef^-^ , (p) which is solved by matrix inversion. 
Then the self-energy reads: 

E„(k,£;) = i^E^'C'i -k)G(ki,S)i^„(ki,k,^). 



(B6) 



The equations (B4,B5,B6) have to be solved self- 
consistently thus determining the self-energy 1]'°' (k, E) . 

Our single-site-approximations neglects all contribu- 
tions from crossed diagrams (leading to weak-localization 
effects, as considered in Ref. |^) as well as the spin- 
resolved electron-electron interaction leading to the split- 
ting of the impurity bands (the Mott transition, see e.g. 
pof). The latter may become important for very low 
densities when the impurity bands are narrow. Within 
this approximation for 17'^*{]s.,E) the integral (^ is a 
well defined quantity, as lm{Y,'^°^ (k, E)} = (and thus 
A(k,^)=0)for£;<i;,„in. 

Finally note that no impurity bands are found within 
the self-consistent Born approximation 



S„(k,i?) - -lE^""(k-ki)G(ki,i?)T/„-(ki -k) 

ki 

(B7) 

which is just the first diagram from Fig. ^ 

APPENDIX C: THE SCREENING 

In order to consider screening we have to include the 
electron-electron interaction given by the Hamiltonian 

^"" = ^11 Wip)aHk + p)aHk'~p)a{k')a{k). (CI) 

k,k',p 

where the matrix element is calculated via 

W{p) = / dzi / dz2**(zi)«'*(z2) 

^e^eop J J 

x-^{z2)'f{zi)e-P^''-'^^ . (C2) 
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Within the random-phase approximation (RPA) the 
screening of the impurity potentials is described by |25 



C(P) 



K(P) 



1 -no(p,w = o)w{p) 



(C3) 



For a free-electron gas the two-dimensional vacuum po- 
larizabihty n"(p, w = 0) for T = is given by ^ 



n°(p,cj = 0) = -po l-e(p-2fcF)Wl-4- 



(C4) 



where kp = {2TrNi:,y^^ is the Fermi wave-vector. 

Actually, the electronic states are affected by the im- 
purity scattering, which may change the density of states 
dramatically as can be seen from Fig. ^. Now the polariz- 
ability n(p = 0) is related to the actual density of states 
at the chemical potential which is significantly lower 
than pq. Calculations within the Born-approximation 
show that the p-dependence of the polarizability becomes 
weaker and that 11(0) decreases with increasing scattering 
In order to accommodate these trends we make 
the replacement n°(fc) n*(fc) = ~p{E^), given by the 
calculated density of states of Fig. || and the chemical po- 
tential at T = 0. Then we obtain the screened impurity 
interaction 



K,(p) 



l + p{EP)W{p) 



(C5) 



This is equivalent to the Thomas-Fermi approximation 
for the screening. The same type of screening has been 
considered in Ref. as well. Of course both ways of 
including screening are approximations. In a full calcu- 
lation the scattering has to be treated self-consistently in 
the calculation of the polarizability. Such a calculation 
was performed in Ref. 1 4^ for a quantum wire within the 
restriction of a delta-potential for impurity scattering. 

Eq. ( |C3| ) only considers screening within the same well. 
The extension to screening by electrons from neighbor- 
ing wells is given in section 6 of Ref. 0|. The results 
are almost indistinguishable for the samples with thick 
barrier width considered in sections IV and V (see also 
Ref. for the screened matrix elements). Screening by 
electrons from neighboring wells becomes more impor- 
tant for a smaller barrier width as used in section VI and 
VII where th e fo rmalism from Ref. |7| was used with the 
polarization (|c|) for RPA and n°(fc) = -p{E^) for the 
TF case. The temperature dependence of the screening 
is neglected in all calculations. 
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FIG. 2. Galculated density of states in units of the 2D free 
carrier density po using RPA screening. The vertical lines in- 
dicate the position of the chemical potential for T = at the 
respective doping densities. 
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FIG. 1. The self-consistent single-site approximation. The 
dashed lines indicate impurity potentials and the double lines 
denote the full Green-function 



FIG. 3. Calculated spectral functions A{E, k) versus 
using RPA screening for different doping den- 
sities at S = 5 meV (a) and E = —5 meV (b). 
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FIG. 4. Calculated temperature dependence of the cur- 
rent-field relations for different doping densities. The screen- 
ing is treated within the RPA. 
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FIG. 6. Calculated density of states in units of the 2D free 
carrier density po using Thomas-Fermi screening for different 
doping densities. 



FIG. 5. Explanation of the two different current max- 
ima within a sketch of the conduction band profile: (a) For 
low temperatures the electrons occupy the impurity band 
(black area). As these states exhibit a flat spectral function 
(see Fig. ^b)) they contain contributions from essentially all 
k-vectors and thus tunneling into the free-particle states is 
possible at all energies. Maximal current is found when all 
states from the impurity band can tunnel into the free-particle 
states, i.e., eFd « |-Emin|. (b) For high temperatures the elec- 
trons occupy the free-electron states as well (the grey scale 
indicates the occupation given by the Fermi-function). The 
spectral function j4(k, E) of such a free-electron state with 
given wave vector k is peaked around E = Ek as shown in 
the figure. Due to k- conservation tunneling can only take 
place if the spectral functions for the same k of both wells 
overlap. On the other hand a net current is caused by the dif- 
ference in occupation. This competition results in a current 
maximum for eFd ~ P as shown in section III. 
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FIG. 7. Calculated temperature dependence of the cur- 
rent-field relations for different doping densities using TF 

screening. 



13 




1 — I — I — r~i — I — I — I — I — I — I — I — I — I — I — I 



-5 
E [meV] 



FIG. 8. Calculated density of states for the sample param- 
eters of Ref. using RPA screening (full line) and TF 
screening (dashed line). 
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FIG. 10. Experimental current-voltage characteristic with- 
out external irradiation together with calculations for differ- 
ent electron temperatures. In the calculation we estimate the 
bias by NFd assuming a homogeneous field distribution and 
neglecting potential drops in the contact regions. 
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FIG. 9. Calculated temperature dependence of the cur- 
rent-field relations for the sample of Ref. ||l^,^. The 
screening is treated within the RPA (a) and within the 
Thomas- Fermi approximation (b). 
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FIG. 11. Temperature dependence of the zero-bias conduc- 
tance for the sample of Refs. [^8|jl^. Full line: Calculation us- 
ing spectral functions from the single-site approximation and 
TF screening, dashed line: Calculation using spectral func- 
tions from the self-consistent Born-approximation, crosses: 
experimental data 
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FIG. 12. Current-voltage characteristics under irradiation. 

a) Theoretical results for hu — 6.3 mcV and different field 
strength eFacd — ahu of the irradiation, b) Experimental 
results for hv = 6.3 meV and different laser intensities in- 
creasing from the top to the bottom. The actual values Fac 
inside the sample are not accessible, c) Theoretical results for 
a — 2.4 and different photon energies. The thin lino depicts 
hv = 5.3 meV and a = 2.1. d) Experimental results for dif- 
ferent photon energies. The laser intensity was tuned to give 
maximum negative conductance. 
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